APS/123-QED 



Takagi-Taupin Description of X-ray Dynamical Diffraction from Diffractive Optics 

with Large Numerical Aperture 

Hanfei Yan, 1, 2, ^j Jorg Maser, 1,2 Albert Macrander, 2 Qun Shen, 2 
Stefan Vogt, 2 Brian Stephenson, 3 ' 1 and Hyon Choi Kang 3,1,4 

1 Center for Nanoscale Materials, Argonne National Laboratory, Argonne, IL 60439 
2 X-ray Science Division, Argonne National Laboratory, Argonne, IL 60439 
3 Materials Science Division, Argonne National Laboratory, Argonne, IL 60439 
^ Advanced Photonics Research Institute, Gwangju Institute of Science and Technology, Gwangju, 500-712, Republic of Korea 

(Dated: February 5, 2008) 

We present a formalism of x-ray dynamical diffraction from volume diffractive optics with large 
numerical aperture and high aspect ratio, in an analogy to the Takagi-Taupin equations for strained 
single crystals. We derive a set of basic equations for dynamical diffraction from volume diffractive 
optics, which enable us to study the focusing property of these optics with various grating profiles. 
We study volume diffractive optics that satisfy the Bragg condition to various degrees, namely flat, 
tilted and wedged geometries, and derive the curved geometries required for ultimate focusing. We 
show that the curved geometries satisfy the Bragg condition everywhere and phase requirement for 
point focusing, and effectively focus hard x-rays to a scale close to the wavelength. 

PACS numbers: 41.50.+h, 07.85. Qe, 61.10.Nz 

I. INTRODUCTION 

X-ray microscopy has found numerous applications in materials sciences, medicine, biology, environmental sciences 
and many other fields by utilizing analytical techniques such as diffraction, x-ray fluorescence and spectroscopy, as 
well as imaging techniques using absorption and phase contrast. Progress in the fabrication of x-ray optics has pushed 
the sp atial re solution of hard x-ray microscopes and -microprobes to well below 100 nm using Kirk Patrick-Baez (K-B) 
mirrors^ 2 * 3 * 4 *, refractive lenses^ and diffractive optics^'. This makes hard x-ray microscopy and spectromicroscopy 
powerful tools to probe the structure and properties of materials at the nanoscale, with the advantages of good 
penetration, which provides for the study of thick specimens and buried structures. As well, there are only moderate 
or no vacuum requirements, enabling in-situ studies, and one can avoid interaction with electric and magnetic fields, 
enabling good quantification. In order to push the spatial resolution of x-ray microscopy to the theoretical limit, 
it is important to know whether there exists a limit in effectively focusing x-rays using optics made of realistic 
materials. Such limits have been studied theoretically for various optics, and yielded numbers of the resolution 
limit of approximately 10 nm for waveguides^ and 2 nm or less for a refractive lens in an idealized adiabatically 
focusing geometry^. For diffractive optics, for example a zone plate with 'flat zones', a recent study found that when 
illuminated by a plane wave a zone plate could have a similar focusing limit as that for a waveguide, while for 1:1 
imaging (spherical wave illumination), where a Bragg condition is satisfied, no limit to the spatial resolution was 
founds Diffractive optics are intrinsically well suited to achieve a high spatial resolution because a large numerical 
aperture (NA) can be achieved by diffraction. A novel approach to making diffractive optics with high NA is the 
Multilaycr-Laue-Lens (MLL) 11 . MLL's are fabricated by multilayer deposition onto a flat substrate and used in Laue 
geometry. An MLL can be considered as a special type of zone plate. The fabrication approach provides high aspect 
ratios, and is therefore well suited for hard x-ray focusing^. Initial 1-D MLL structures have achieved a line focus of 
30 nm in the hard x-ray range (19.5 keV) with a diffraction efficiency of more than 40% 6 . 

MLL's exhibit strong dynamic diffraction properties, and have to be described using a dynamical diffraction theory. 
Using such an approach, namely Coupled- Wave-Theory (CWT^Pl, we have shown that MLL's can reach a spatial 
resolution of one nanometer^. However, these simulation results were limited by assumptions made in the derivation 
of the local one-dimensional (1-D) CWT equations for volume diffractive optics^. To understand the ultimate limit 
of focusing by diffractive x-ray optics made of realistic materials, we have developed a new theoretical approach 
to volume diffraction, based on the Takagi-Taupin description of dynamical x-ray diffraction developed for strained 
single crystals. This approach represents a full- wave dynamical theory for volume diffractive optics with high numerical 
aperture, ft takes into account flat, tilted, wedged and curved zones as well as arbitrary zone profiles, allowing us 
to simulate any of these geometries for imaging and focusing. Since it is based on a dynamical diffraction theory, 
it provides physical insight into the diffraction properties of volume diffractive optics. Using this approach, we have 
studied the focusing performance of various MLL's, and have investigated limitations that constrain the achievable 
NA. We find that a wedged MLL is able to focus hard x-rays to less tahn a nanometer, while curved zones are needed 
to achieve a beam size approaching the x-ray wavelength. 
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II. THEORETICAL APPROACH 

In the conventional geometrical-optical theory, a zone plate is considered to be "thin" so that the volume diffraction 
effect can be neglected^. This approximation is only valid when the zone plate thickness, w, is less than (2Ar m i n ) 2 /X, 
where Ar rn i n is the outmost zone width and A is the wave length of the incident x-ray. Such a condition is usually 
satisfied for diffractive optics at optimum thickness with a zone width no smaller than 10 mrP^. A more rigorous full- 
wave theory that takes into account the dynamical diffraction effect is needed to correctly describe x-ray diffraction 
from a zone plate with outmost zone width less than 10 nm because the optimized thicknesses are too large to treat 
them as "thin". This was initially done by Maser and Schmahl 14 , who applied locally a 1-D approach of CWT to 
study the variation of the local diffraction efficiency with thickness and with the slant of the zones. Their approach 
assumes that the zone plate can be decomposed locally into periodic gratings with fixed gratings period. Their 
approach properly accounts for dynamical diffraction effects, and allows the calculation of the point spread function 
and modulation transfer function of zone plates with high aspect ratios and small outermost zone width. However, 
their approach is limited to cases of w << / and relatively small NA, where / is the focal length of the zone plate. 
For an optimized zone plate thickness (these are less than 1 /im for soft x-rays and greater than 1 /im for hard x-rays) , 
and for a zone plate radius matched to the coherence length of a 3 rd generation synchrotron source, their approach is 
limited to MLL's with outmost zone width no smaller than 1 nm. 

Other calculations have been published since that of Maser and Schmahl. The 1-D CWT app roach was further 
extended to study the effects of an arbitrary line/space ratio and interface roughness by Schneide j 16 * 17 l By decom- 
posing zone plates into local gratings, and with similar limitations as 1-D CWT, Levashov and Vinogradov^! studied 
the variation of the total diffraction efficiency with thickness. A different numerical approach based on the parabolic 
wave equation within paraxial approximation has been employed by Kurokhtin and Popov^H Their numerical method 
requires large computing power when thousands of individual zones ar e con sidered. More recently methods of solving 
the parabolic wave equation using eigenfunctions have been reporte d 10 * 20 !. Due to the use of the paraxial approxi- 
mation, numerical approaches based on the parabolic wave equation are limited to intermediate values of numerical 
aperture (NA). To overcome the limitations of these theories, and to provide an approach that is valid to wavelength 
scale, we have developed and present below a new approach to modeling dynamical diffraction from MLL's based on 
first principles. 

From Maxwell's equations, we obtain the scalar-wave equation describing the electric field variation of a monochro- 
matic x-ray wave in a medium with susceptibility xif) as i 

V 2 E{r) + k 2 [l + X {r)]E{r) = Q (1) 

where k — 2ir/\ and E is the electric field vector. In principle, for a given function x(r), the electric field E can be 
solved directly from Eq. (JT|). However, in many cases the variation of x(f) 1S t°° complicated to allow an analytical 
solution to Eq. (1). Even using a numerical approach, the computing time can be prohibitive. An extreme case 
is x-ray diffraction from single crystals, where tens of thousands of atomic planes are involved in the scattering 
process and solving Eq. ([IJ plane by plane becomes impractical. The same difficulty arises for x-ray diffraction 
from volume diffractive optics, which contains thousands of diffractive structures, and requires the consideration of 
the multiwave scattering process. In crystallography, this difficulty has been overcome elegantly by the classical 
dynamical diffraction theory, in which the periodic susceptibility function x is expanded into a Fourier series. Each 
Fourier expansion coefficient, Xh, serves as a resonator, interacting with the incident plane wave and exciting a 
resonant wave. The excited resonant wave will again interact with resonators and excite high order resonant waves, 
leading to a dynamical multiwave scattering process and a self-consistent system of equations. Instead of seeking a 
solution to the second-order partial differential equation [Eq. ([!])] for each atomic plane, the diffraction problem in a 
single crystal is simplified to finding a set of plane waves that satisfy the self-consistent system, thus allow a much 
simplified mathematical and numerical solution. To deal with strained single crystals, this approach was generalized 
by Takagi and TaupirP32U. An additional phase term is added to each Fourier component of the susceptibility function 
to account for the phase change caused by the deformation of the lattice planes due to misorientation and strain fields. 
Detailed treatments of x-ray dynamical diffraction theory can be found in many monographs^ 2 ^. We show that a 
similar methodology can be applied to solve the diffraction problem for diffractive optics. 

We begin by considering a Fresncl Zone Plate (FZP). A FZP is a circular diffraction grating that is capable of 
focusing an incident wave to a series of focal points. It consists of alternating zones made of two different materials, 
and, in the limit of a 'thin' zone plate, produces a phase shift tt between waves passing through neighboring zones, 
thereby yielding a focusing effect. The position of the nth zone is determined by the zone plate law, 



r 2 n = n \f + n 2 A 2 /4, 



(2) 
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where the width of the nth zone is given by, 




(3) 



It is seen from Eq. d3| that the zone width Ar„ decreases radially. Locally, the zone width does not vary rapidly. 
When r n is large and^r„ is small, a FZP might be viewed as a "strained crystal" with c£-spacing equal to d = 2Ar n . 
Therefore, we reason that a formalism akin to the Takagi-Taupin equations in crystallography can serve to allow 
a study of dynamical diffraction effects in a volume zone plate or MLL, and, more generally, in volume diffractive 
optics. A difficulty arises when applying the Takagi-Taupin equations since the d-spacing of volume diffractive optics 
changes significantly from its center to its outmost region. Its susceptibility function can therefore not be expanded 
into a Fourier series in a conventional way. We address this by finding a mapping relationship between zones in a FZP 
and zones in a strictly periodic structure with a known Fourier series expansion. From this mapping relationship, we 
expand the susceptibility function for a FZP into a pseudo Fourier series. 

We will study in this paper flat, tilted, wedged and curved geometries (see Fig. [I]). For the sake of simplicity, 
we begin by considering a 1-D MLL with flat zones. These zones can be aligned parallel to the optical axis in 'flat 
geometry', satisfying the Bragg condition on the optical axis, or in tilted geometry, satisfying the Bragg condition in 
some outer area of the MLL. The zones are made of materials A and B, whose layer thickness follows the zone plate 
law Eq. ([2]). This geometry has been manufactured as MLL's with outermost zone widths of 15 nm, 10 nm and 5 nm, 
and demonstrated to focus hard x-rays efficiently into a line with a width smaller than 20 nm 25 . The aspect ratios 
of the manufactured MLL structures range from 1000 - 3000, more than an order of magnitude larger than what 
has been achieved using lithographic techniques. To find the pseudo Fourier series of an MLL, we consider another 
periodic binary structure with mirror symmetry about the axis z as shown in Fig. [2j We consider that it is made 
of the materials A and B and has a A/B layer ratio of 1:1. This periodic grating has a radial coordinate, x', that is 
different from the radial coordinate, x, for the MLL. A connection between x and x' will be established later. The 
Fourier series of a periodic structure can be written straightforwardly as, 

XP (x') = * + f^[l- (-1)^(2/^) 



E 



A X n , l\|h|l fni._\ X ' 



, 2ihn 

h— — oc,h^0 



[l-(-l)W]exp(2/ l7 rV): ^0 (4) 



where \ — (XA + Xb)/2 , Ax = Xa ~ Xb and T is the period. Here T is taken as arbitrary. If a relationship between 
x and x' can be found so that the nth layer in the perfect periodic structure is mapped to the nth layer in a 1-D 
MLL, a series representation of the susceptibility function for the MLL can be obtained as, Xmll(x) = Xp[ x '( x )}- 
From the zone plate law [Eq. Q], the nth layer index, n, can be written as, 

n=\(JxlVp-f). 

Similarly, in the periodic structure the nth layer index, n', can be written as, 

/ = 2_KJ 
T 

By equating n and n' we establish a relationship between x and x' given by, 

W\ = j(V xY Tp-.f). 

Substituting this into Eq. Q we obtain the following series expansion, 

00 

Xmll(x)=Xq+ E Xhexp[ic/) h (x)}, x^0, (5) 

h=—oo,hjt0 

where xo = Xi Xh = (Ax/2i/wr)[l — (— l)' fc '], and fa = hk(\j 1 x 2 + f 2 — /). Since the Fourier series of the periodic 
structure is mathematically complete and x' is one-to-one mapped to x, the series expansion in Eq. §5§ can be 
considered as a pseudo Fourier series. Using Eq. (pj, we can construct the MLL structure, aside from the singular 
point at origin. For simplicity in the following discussion the subscript notation of MLL of xmll will be dropped. 
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We note that the gradient of the phase function, Vifih, is equal to the local reciprocal lattice vector, ph = (2irh/d)e x , 
where e x is a unit vector along the positive x axis. Unlike the conventional Fourier series representation of a strictly 
periodic structure, here the zone plate structure is represented by a superposition of sine waves with varying periods 
obeying the zone plate law (see Fig. [3]). If the radius of an MLL is sufficiently small so that x max is less than 
(4/ 3 A/ l^-l) 1 / 4 , the following approximation can be made, 

x n ~ nXf, and <j> h ~^x 2 . (6) 

The basic zone plate law is thereby reduced to the first term in Eq. Q only. 

For an incident wave with an arbitrary wave front, we assume that it can be written as, 

Einddent = Eq if) exp(lfcs ' f), (7) 

where sq is a unit vector along the incident propagation direction, and E^ is a slowly varying complex function that 
represents the modulated amplitude of the incident wave. When an MLL is illuminated by such an incident wave, we 
write a trial solution to the wave equation as, 

E = 4(r) exp^fcso • r + <j> h )) = £ E h P h , (8) 

h h 

where Ph — exp[i(kso ■ f+ 4>h)]- This scries form of the trial solution is chosen in accordance with the zone plate 
structure and Eq. ^ . Because a zone plate structure can be represented by a pseudo Fourier series, the phase term 
Ph in the trial solution in Eq. (JsJ) indeed reflects the primary component of the spatial variation of the wavefield inside 
the global structure. Recall that ph — V^/,, so in a local region the trial solution in Eq. Q is a Bloch wave solution. 
In addition, a further examination of Ph shows that it actually corresponds to the phase variation of a spherical wave 
that converges to the hth order focus of the MLL, therefore this trial solution represents a superposition of a set of 
spherical waves with modulated amplitudes that converge to MLL's foci. 
Substituting Eqs. ^ and Q into Eq. Q we obtain, 

[^E M -VPh + PhV 2 E hl + E hl (k 2 P h + V 2 P h )]e t ^-k 2 J2Ph(J2^iEi), 

h i—x.y,z h I 

where e x .y.z are unit vectors along x, y and z axis. By equating them for each we obtain an infinite set of equations, 
J2 [2VE M ■ VP h + P h V 2 E hl + E hl (k 2 P h + V 2 P h )]e t = ~k 2 P h Xh-lEi, (9) 

i—x,y,z I 

M = 0,±l,±2,±3,... 

If there is a solution to Eq. (10 1, it will be also a solution to Eq. (JT|). Eq. (10 1 can be simplified using the following 
assumptions. First, the second order derivative, V 2 Ehi, is negligible compared to other terms. This is equivalent to 
saying that the amplitude envelope function, Eh , varies very slowly over a length scale equal to the wavelength or 
zone width, i.e., only the phase term, Ph, contains quickly varying components. For this approximation to be valid, 
one needs to ensure that no specularly reflected wave from the MLL/vacuum boundary is strongly excited. This is 
because if the incidence angle (the angle between the incidence direction and the MLL surface) , is much larger than the 
critical angle at which total external reflections occurs and the specularly reflected wave is negligible, the approximate 
solution obtained without second order derivatives can be uniquely determined by the boundary condition for the 
electric field alone. Under this condition, it is easy to verify that to a good approximation this solution can also 
satisfy the boundary condition for the magnetic field. In other words, without second order derivative terms, the 
system is still self-consistent. On the other hand, when the incidence angle approaches the critical angle and the 
specularly reflected wave is strongly excited, this approximate solution can no longer satisfy the boundary condition 
for the magnetic field, violating the physical law for electromagnetic field. To have a solution that is able to satisfy the 
boundary conditions for both electric and magnetic fields, second order derivatives should be retained in the equation. 
For hard x-rays the critical angle for most materials is smaller than 1°. For example, at 20 keV the critical angle of 
total external reflection is 0.23° for tungsten. The specular wave decreases below 10 -3 above 0.67°. If we consider 
values of the numerical aperture of 0.7 or smaller here, corresponding to incidence angles of 45° or larger, then the 
second order derivatives can be neglected. It is important to point out that this approximation of neglecting second 
order derivatives also impose a restriction on how rapidly E^ 1 can vary on the entrance surface, in order to guarantee 
that the second order derivative terms are negligible everywhere. 
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By neglecting second order derivatives we obtain, 

Y (2V£ W ■ VP,, + E hl X7 2 P h + k 2 E hi P h )ei = -k 2 P h Y Xh-lA 

i—x,y,z I 

Taking the dot product on both sides with one obtains a scalar equation, 

VE 2 h ■ X7P h + E 2 X7 2 P h + k 2 E 2 P h = -k 2 P h Y Xh-iE h Ei costf w , 



(10) 



where dhi is the angle between the polarization directions of two wave components, Eh and Ei. Eq. (10 1 can be 
further simplified by canceling one factor of E^ to give, 



2VE h ■ \7P h + E h V 2 P h + k 2 E h P h = -k 2 P h Y Xh-lEt cos W . 



(11) 



Because P h is known [see Eq. d8j], we obtain, 



VP ft = i(k + V(f> h )P h , k = ks , 



V 2 Ph 



iV 2 ^ - [k 2 + 2k • V<t> h + {V<f> h f]P h . 



Substituting these expressions into Eq. (Ill yields 
2i 



VE h ■ (s + ^) + f3 h (r)E h + Y Xh-lEi cos W = 0, h, I = 0, ±1, ±2, ±3, . . . , 



fc 



(12) 



where 



k 2 



2s ° nr 



k ' 



k 2 - (fcp + Vcj> h ) 2 
k 2 



These are the central equations that will be using in this paper. Here f3h is the deviation function, which quantifies the 
violation of the Bragg condition. Because an infinite number of differential equations need to be solved simultaneously 
(both h and I run from -oo to oo), a system described by Eq. (12) is considered presently to be still too complicated. 



By noting that a diffracted wave is strongly excited only near the Bragg condition, we can truncate the system to a 
finite number of equations. Consequently, it is sufficient to consider a finite number (N) of coupled first order partial 
differential equations, with the boundary conditions of the electric field at the entrance surface (Laue case)given by, 



Eh,hjto — 0, 



E = E ( a) - 



(13) 



One may recognize that Eq. ( 12 ) resembles the Takagi-Taupin equations for strained single crystals. This is due to 



the fact that the methodologies of these two theories are the same: both utilize the prior knowledge about the solution 
for diffraction from a periodic structure, thereby allowing the separation of the fast and slowly varying components. 
However, there is a fundamental difference between them: the lattice constant of a strained single crystal differs 
from its unstrained value only very slightly, and as a result the diffracted wave vector can be regarded approximately 
invariant with position. In other words, in a strained crystal, the following approximation is valid: 



so 



so 



Ph 
k 



so 



Ph 

k 



Sh, 



(14) 



where ph is the unstrained reciprocal lattice vector and Sh is a unit vector along the diffracted wave direction. In 
addition, in a crystal often only one diffracted wave is strongly excited, so only two waves need to be considered. 
Under these approximations, Eq. ( 12 1 becomes, 



£- o E =il(xoE + Cx h E h ) 
£- h E h = if [{p h + X0 )E h + CxhEo] 



where C = cos d 



1)0: 



(15) 



which are the well-known Takagi-Taupin equations for crystal diffraction. For volume diffractive optics its d-spacing 
can vary from hundreds of nanometers in the vicinity of the center to several nanometers or even smaller in the outmost 
region, so the variation of the diffracted wave vector is significant and the approximation in Eq. (14.1 is invalid. Also, 



for areas of a diffractive optics with large d-spacing, many diffraction orders can be excited, and consideration of only 
two orders is insufficient for the description of volume diffractive optics in these areas. 
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III. DIFFRACTION FROM AN MLL WITH FLAT, TILTED AND WEDGED ZONES 

A. Flat and Tilted Zones 

Using a Takagi-Taupin description for volume diffractive optics, we can study dynamical diffraction from various 
types of MLL's. We will discuss flat, tilted and wedged MLL's. Flat and tilted MLL's have zones that are parallel to 
each other, but are aligned at different tilt angle with respect to the incident beam. Wedged MLL's have zones that 
are tilted with respect to each other, such as to satisfy the Bragg angle with regard to the direction of the incident 
wave. We consider MLL's with a diameter of 30 /im, which corresponds to typical lateral coherence lengths of 20 
keV x-rays at a 3 rd generation synchrotron source, and a thickness of 13.5 /im, which corresponds to the optimum 
thickness of an MLL of Si and WSi2 for a photon energy of 19.5 keV. As an example, we use an outermost zone width 
of 5 nm, corresponding to MLL structures we fabricated currently^. A plane wave at 19.5 keV is impinging on the 
MLL described above, with an inclination angle of 9 to the normal direction of the sample surface, as shown in Fig. [4] 
For simplicity, cr-polarization, corresponding to cosi?w = 1, is assumed. In this case x max — 30 /im is much smaller 
than the focal length / — 4.72 mm and the approximate zone plate law can be used, so, 

a? 1 ' 



x n = \J n\f, and <p h = 



From Eq. (12 1 one can deduce, 



/3,(*)«-2^sin0-(^) 2 , 
V<^/, hx _ 

In order to obtain a complete description of the wavefield, 11 beams (h = 0, ±1, ±2, ±3, ±4, ±5) are considered. Thus, 
we need to obtain a solution for eleven coupled hyperbolic partial differential equations given by, 

[(sm6+-x)— + cos 6— ]E h = ij(Y, X h-iEi + h E h ), h = -5, -4, ...4, 5, (16) 

■I X ' i = - 5 



with the boundary conditions, Eh^o(x,0) = and Eq(x, 0) = 1. Eq. (16 1 is solved numerically in Mathematica 5.2, 
and the simulation results are displayed in Fig. 5(a)| for 9 equal to 0, 1.6 and 3.2 mrad., respectively. Since we are not 



interested in the diffractions with positive orders which correspond to divergent waves, they are not plotted. We find, 
as shown in the figures, that higher order diffractions have weaker intensities and the fifth order diffraction intensity 
becomes negligible. This justifies the approximation of only retaining diffraction orders up to order five. 

For normal incidence (9 — 0), the intensity distribution of all diffracted beams show a similar trend: only in 
the vicinity of the center (x — 0) do the diffracted beams have appreciable intensities; below a certain d-spacing the 
diffracted intensity decreases to zero quickly. The local diffraction efficiency, defined as the ratio of the local diffraction 
intensity to the incident wave intensity, which is unity, is below 30% for all diffraction orders. Near the center where 
the value of the deviation function, /3/j, is relatively small, a large number of diffracted waves are excited. The value 
of the deviation function increases rapidly as the e?-spacing becomes smaller with increasing radius, causing a fast 
drop of the diffraction intensity in the outmost region. This indicates that only a fraction of the MLL near the center 
contributes to the focusing. This is consistent with the results obtained by Maser^l and by Pfeiffer et al.™ The 
radial variation of the diffraction efficiency also depends on the thickness of the MLL, as seen in Fig. |5(b)| where the 
— I s * order diffraction intensity profile is plotted as a function of radius and thickness. When the thickness is very 
small, the diffraction efficiency is low, and no dynamical effects are encountered. At larger thickness, the efficiency 
increases for areas with larger d-spacing, but stays small for areas with smaller d-spacing, where deviations from 
the Bragg condition have an appreciable effect. This trend is clearly shown in Fig. |5(c) where the radial efficiency 



distribution of the —1 st order at different thickness is plotted. A maximum diffraction efficiency, about 34%, is reached 
at w = 10 /mi. This number agrees with the value of optimum thickness for a phase zone plate obtained in geometrical 
optics, w — A/Ax ~ 10 /LtrrPSI 

When the MLL is tilted, the Bragg condition can be satisfied for a small area of the MLL. The middle graph in 



Fig. 5(a) shows the local diffraction efficiency at a tilt angle of 9 = 1.6 mrad. At this angle the Bragg condition is 
satisfied for the —I s * diffraction order at x — 15 /im, corresponding to a zone width of 10 nm. A diffraction efficiency 
maximum is observed at x = 15 /im, with a local diffraction efficiency of 67%, as expected from dynamical diffraction. 
At this tilting angle the local diffraction efficiency starts to decrease at a larger radius than that in the case of normal 
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incidence, indicating that a larger fraction of the MLL contributes to the focusing. An increase in the diffraction 
efficiency of the —2 nd (focusing) diffraction order is observed for a radius of 7.5 /mi, corresponding to a zone width of 
20 nm, where the Bragg condition is satisfied for the —2 nd order. In Fig. 5(b) the variation of the local efficie ncy of 
the —1 st order diffraction as a function of radius and thickness is displayed. The bottom graph in Fig. 5(c) shows 
its radial efficiency distribution at different thicknesses. We observe at a thickness of 16 /iui, a maximum diffraction 
efficiency of the —1 st order is reached. In the phenomenon of "Pendellsung fringes" for dynamical diffraction from a 
single crystal in Laue geometry, a Laue peak maximum is reached at a thickness of w — A cos 9/2y/\x-iXi\ ~ 16 ^tirPH, 
which is in a reasonable agreement with the value obtained here. This agreement suggests that dynamical resonant 
scatterings are strongly excited in this case, which raises the diffraction efficiency as well as the effective radius. 

Here the even order diffractions are excited via dynamical multiwave scattering processes, similar to the Renningcr 
effect in diffraction from crystals^, although they are forbidden in the geometrical-optical theory^. 

To study the focusing property of the flat and tilted MLL's, we employed Fresnel-Kirchhoff integral to calculate 
the wavefield at a point behind the MLL. In Fig. 6(a) the isophotes near the focus at different tilting angles are plotted. 
As expected, we observe in our calculations that tilting leads to a reduction in beam size and to an enhancement of 
peak intensity. However, there is an optimum tilting angle above which defocusing occurs. This tendency is clearly 
seen in Fig. |6(b)| by comparing the line focus profiles on the best focal plane at different tilting angles. When 8 = 0, 
the full width at half maximum (FWHM) of the best focused beam is about 7.4 nm. As 8 increases to 2.1 mradian, 
the FWHM of the best focused beam is reduced to 4.7 nm because a larger fraction of the MLL contributes to the 



focusing, as shown in Fig. 5(a) However, when 8 = 3.2 mradian at which Bragg condition is satisfied for the outmost 
zones at x — 30 /im, multiple peaks with nearly equal intensities appear on the plane of best focus, destroying the 
resoluti on. T his is a result of the uneven radial intensity distribution at 8 = 3.2 mradian (see the bottom graph 



in Fig. 5(a) I. At this tilting angle, the diffraction efficiency of the —I s * order has appreciable values only near the 



center and near the outmost regions where a Bragg condition is satisfied, while in between the diffraction efficiency is 
very low. As a result, the diffracted waves from the center and the outmost region interfere with each other strongly, 
causing multiple peaks with nearly equal intensity on the focal plane in an analogy to the multi-slit diffraction pattern. 
The effective radius within which the diffraction efficiency is nonzero constrains the focal size of a tilted MLL with 
parameters given here to slightly smaller than 5nm. Therefore, if the radius of the MLL is increased to 60 /jm and 
the outmost zone width is reduced to 2.5 nm, at this optimum thickness we still cannot obtain a 2.5-nm beam size. 
To increase the effective radius, we need to reduce the thickness of the MLL, but with a significant loss of efficiency, 
as shown in Fig. |5(c) 



The improvement of the focusing by tilting might be better understood by considering the structure of an MLL, 
which consists of many thin layers. The diffraction property of an MLL is determined by the scattering property of 
a single layer and the interference effect of the scattered waves from different layers as well. Each thin layer can be 
regarded as a waveguide where total external reflection occurs at the boundaries. As Bergemann et al. arguecP, the 
critical angle of the total external reflection imposes an ultimate limit on the achievable NA for a waveguide, because 
above this critical angle x-rays start to leak out, resulting in a rapid drop of the reflection intensity. Consequently the 
maximum converging angle of x-rays from a waveguide can not exceed the critical angle, so the maximum achievable 
NA for a waveguide is about v / 2An, where An is the difference of the refractive index of the waveguide material from 
unity. They also claimed that this limit applies to all x-ray focusing optics. However, it has been shown that this 
limit could be overcome by a thick adiabatic refractive lens^. Such a limit does not apply to diffractive optics either 
since resonant scattering occurs when the Bragg condition is fulfilled, so that the reflected waves from different layers 
can interfere constructively, leading to a resultant wave that can be as strong as the incident wave. For the MLL the 
converging angle is only limited by the Bragg angle, which can be as large as w/2. Though the waveguide effect will 
somewhat influence the achievable NA of an MLL when Bragg condition is not satisfied (8 = 0), it docs not determine 
the ultimate resolution of an MLL. 



B. Wedged MLL 

As shown previously, tilted MLL's can focus x-rays much more effectively than flat MLL's. However, tilted MLL 
cannot achieve a 1-nm focusing without significantly sacrificing the efficiency. To achieve both high efficiency and 
small focusing, we consider a structure in which there is a distinct tilt to every zone by an angle increasing with radius 
so that the local Bragg condition is fulfilled in every zone. This geometry corresponds to a wedged zone shown in Fig. 
[I] Deducing from Eqs. ^ and (|3j, we obtain the modified zone plate law corresponding to an MLL with wedged 
zones which satisfy the Bragg condition at z = 0, 



Ti A^ z 

a„(z)\lnXf-\ , a n = 1 ; , — r. (17) 

1 J 4 ' 2/(1 + tiA/4/) V ; 
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At the entrance surface (z = 0) the zone width obeys the conventional zone plate law [Eq. p|], but all zones shrink 
by a factor a n (z) along the depth z. The shrinkage factor, a n , also depends on the zone number n. To simplify 
the calculation, as an approximation we consider a shrinkage factor independent of n given by, a n = 1 — z/[2f(l + 
rWcA/8/)], where n max is the maximum zone number. Consequently all zones shrink homogeneously by a factor 
proportional to the depth z, 



and the phase function (p^ becomes, 



4> h = hkU 



a(z)'- 



+ f 2 -f], 



h = 0,±1,±2,±3, 



(18) 



To explore a possible upper limit of the achievable NA, we consider an MLL with following parameters: / = 4.72 
mm, w = 13.5 /jm, E = 19.5 keV, x max — 600 /im and an outmost zone width of 0.25 nm. As we demonstrated 
previously [see Fig. 



5(a) 



for zone widths less than 10 nm high diffraction orders are negligible and the two-beam 
approximation is sufficient to describe the diffraction property of an MLL. Since here we are interested in the diffraction 
behavior of zones with much smaller widths, from now on, unless otherwise specified, the two-beam approximation 

and [7(b)] In Fig. 



will be used. Simulation results are shown in Fig. 7(a 



7(a) the radial efficiency distribution and 



the phase deviation of the diffracted wave from a perfect spherical wave that converges to the focal point are plotted. 
Although in the outmost region the diffracted waves still have strong intensities, they do not contribute to the focus 
because their phase deviation varies rapidly over it and they interfere destructively at the focal spot. Therefore, using 
this structure, we cannot focus an incident plane wave down to the diffraction limit due to the phase effect. In Fig. 
|7(b)| the intensity profile on the best focal plane is displayed. The inset on top shows the isophotes near the focus, 
and the bottom one is a sketch of the zone plate structure. In this case the achievable NA is also limited by the 
dynamical diffraction property of the MLL and a larger physical radius will not increase its NA. This study indicates 
that to effectively focus x-rays down to the wavelength using diffractive optics, not only the Bragg condition has to be 
satisfied so that the diffracted wave intensity is strong and the efficient is high, but also the phase has to be right so 
that all diffracted waves add up in phase at the focal spot. We will derive the optimum curved shapes that minimize 
the radial phase change in section IV. 



C. Localized One-Dimensional Theory 



In this section we show that under certain conditions the assumption of local perfect gratings becomes valid and 
the basic diffraction equations derived in this paper are simplified to simpler forms similar to those in 1-D CWT. For 



the sake of simplicity, we start with the two-beam approximation, so that Eq. (16 1 is reduced to, 

/i = ±1,±3, ±5, . . 



(sin0£ +cos0£)£ o = if (x E + x- h E h ) 
[(sin0 + hx ^ d 



T'di -TCosV£]E h = if[(xo + Ph)Eh + XhEo 
To simplify the calculation, dimensionless variables are used, 

hx hz 



(19) 



/sin0' 



/cos I 



Consequently we obtain, 



[(1 + *)lk + M E h = ^[(Xo + Ph)E h + XhE ] ' 



(20) 



where f3h ~ — x(2 + x) sin 0. We may further simplify Eq. (20 1 by setting 



~ f 77 
E = E exp(i—xoz), 

~ f 77 
E h = E h exp{i—xoz), 
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so the refraction and photoelectric absorption effects are taken out of the equation. Substituting above expressions 



into Eq. (20 1 yields 



[(! + i + II^ = i&(0hEh + XhE 

which are two coupled hyperbolic partial differential equations. According to the partial differential equation theory^, 
a solution at a point (io, Zq) to a hyperbolic equation only depends on the initial data in a bounded domain of 
dependence; the outside part of this domain has no influence on the value of the solution at (xq,zq). For a system 



described by Eq. (21 ), the domain of dependence is bounded by its characteristics, 

x = (x + 1) exp(z — z ) — 1, x = x + (z — zq). (22) 

In Fig. [8] the domain of dependence of point P is depicted. The counterpart concept in crystallography is the inverse 
Borrmann triangle in a single crystal, where the wavefield at P is determined by the wavefields in the triangle RPQ 
(Fig. [8]). A physical elucidation of this concept is that because waves propagate along either so or sh, the wavefield 
at P can only be influenced by the wavefields at points that can propagate to P. Such points form a domain of 
dependence bounded by RPQ. In single crystals both sq and Sh are invariant, so the shape of the inverse Borrmann 
triangle is independent of position. However in an MLL the diffracted beam direction Sh is a function of position 
because the change of reciprocal lattice vector causes a curved edge QP which will vary with position. 
By noting that for w << /, z << 1, one can make the following approximation, 

exp(z — Zq) W 1 + Z — Zq. 



Then the first equation in Eq. (22 1 becomes, 

X = Xq + (x + 1)(Z - Zq). 



The physical manifestation of this approximation is to replace the curved edge QP by a straight line, so Eq. (21 ) can 
be simplified. With the replacement of variables, 

£ = z — x, 

T) = X — (x + l)z, 

one obtains, 

d d d 

dx dz dr/ ' 

_ d d d d 

(1 + x)— + — = -x— - (1 + x)z — . 
ox oz at; orj 

In case is very small so that z « x can be satisfied, the following approximations are valid, 

— — - -x— 
dx dz dr] ' 

. d d d 

{l + x) dx + dz K ~ X W 



Use of these approximations leads to, 



(23) 



If we note the fact that the wavefield at P is only affected by the wavefields inside the triangle RPQ, and under 
conditions z << x and z << 1 the value of x and (3h change only very slightly inside this triangle, it is justified to 
replace variables x and (3h by their values at P. In this way we arrive at two equations that can be solved analytically: 



d_ 
/'// 

d 



Eq = iXh E h(€,v) 



^E h = i/3 h ($, Tj)E h (£, rj) + iXhEo(Z, f?) ' 



(24) 
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where, 

~ _ _ /?r Xh.h 
Xh ' h ~ hX x(P) ' 

~ fn[k(P) 
l h hX x{P) ' 



are constants. One may recognize that Eq. (24 1 coincides with Takagi-Taupin equations for perfect single crystals, 
implying our approximations have the following physical meaning: locally an MLL can be treated as a strictly periodic 
grating. But we want to emphasize that the validity of these approximations is limited to the case where is small 
(corresponding to a small NA), and w « f, i.e., within the triangle RPQ the variation of periodicity can be 
neglected and MLL can be considered as a periodic grating locally. We found that there were similar discussions 
about the validity of local 1-D CWT in volume holography previousl} * 30 * 31 [ If all beams are retained in the equation, 



the discussion is the same but Eq. (24 1 will become N-beam Takagi-Taupin equations for perfect single crystals. From 



Eq. (24 1, a close- form analytical solution can be obtained. 



IV. IDEAL STRUCTURES FOR HIGH NUMERICAL APERTURE 

We have shown above that a wedged MLL cannot achieve a focus close to the wavelength. Ideal structures, which 
not only satisfy Bragg condition everywhere so that the diffraction efficiency is high, but also have a right phase 
locally so that all diffracted waves interfere constructively at the focus, are needed. The first condition requires the 
deviation function, to be zero everywhere. We may consider this condition in reciprocal space with the aid of 
Ewald sphere. In Fig. 9(a)| we show a case in which a plane wave is incident normally on a zone plate with local 



reciprocal lattice vector of ph- For a Bragg condition to be satisfied at a point (x, z) within the MLL, the end point 
of the vector kh — ko + ph should lie on the Ewald sphere, which has a radius k and is centered at (x,z), so that 
j3h = 0. In addition, for the diffracted wave to converge to the focus F, all kh should point to it. In order to satisfy 
both conditions, kh has to take the form given by, 

r ~kx _ k(f — z) _ 

Vx 2 + (/ - z y v^ 2 + (7 - 

Because ko = ke z , ph can be determined unambiguously from the relationship kh — ko + ph- Then, according to the 
relationship, ph — V4>h , we can obtain <ph by integrating ph- We then obtain the following equation, 

4>h = — k[-\/x 2 + (/ — z) 2 + z + const] . 

The integration constant in (ph can be arbitrary, because it does not affect the Bragg condition. This constant is 
determined by the second condition: all diffracted waves should interfere constructively at the focus F. By calculating 
the optical path, we obtain, 



b h = -k[y/x> + U-z) 2 - (/ - z)\. (25) 



We can see that only the —I s * order diffraction (h = —1) can satisfy both conditions. From Eq. (25 1 the ideal zone 
plate law for an incident plane wave can be deduced, 



- 2 - nX(f - z) + n 2 X 2 /4, (26) 



which is a set of confocal paraboloids [Fig. |9(b)| . Substituting Eq. ( [26] ) into Eq. (12), with the two-beam approxi- 
mation we can simplify the system to two partial differential equations, 



d -E = if(xoE + XiE- 1 ) 



{ ^+(/-,)' +fr - *) & E ~i = ^(X-iEo + Xoig-i) • (27) 

Solved numerically, the efficiency and the phase deviation of the —I s * order diffraction that we obtain are shown in 



Fig. 10(a) High diffraction efficiency, 67%, is observed over all radial positions, and the phase deviation is rather 
small, indicating that a converging spherical wave front with very small distortion forms on the exit surface of such an 
MLL. The abrupt drop of the diffraction intensity close to the outmost boundary is a result of incomplete structure 
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beyond it. In Fig. 10(b) we show the intensity profile on the best focal plane with a FWHM of 0.21 nirPH, anc j 



the inset on top shows the isophotes near the focus. Comparing to the peak half width of 0.34 nm obtained from a 
wedged MLL with a same outmost zone width of 0.25 nm, we conclude that curved zone profiles are needed to achieve 
a resolution approaching the wavelength. 

In a similar way, the ideal structure to focus a spherical wave (from a point source) can be obtained as well [see 
Fig. 11(a)] . The difference is that in this case the incident wave vector is a function of position too, 



k = 



k.i 



k(l + z) 



^X 2 + (T + z) 2 y/x 2 + (l + Zf 



where l is the distance from the point source to the entrance surface of the zone plate. To focus it to a point located 
at a distance U away from the entrance surface on the downstream side, we find, 



<p h = -k[y/x 2 + (T + Z) 2 + y/x 2 + (Zj - Z) 2 - (l + h)\ , 



then we obtain the zone plate law for focusing a spherical wave, 

Ax 2 n A[z+(l Q 



n 2 \ 2 /A + n\(l + U) (n\/2 + l + k) 2 



(28) 



(29) 



Again, only for h = — 1 can all requirements be satisfied. Eq. p9|) describes a set of confocal ellipsoids, as shown in 



Fig. 11(b) 



It is interesting to consider 1:1 imaging in this case. If l + h = 4/ and l Q — h = —w, Eq. (29 1 then becomes 



nXf- 



n 2 X 2 n 2 X 2 /A + AnXf 
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(nA/2 + 4/)^ 



{z-w/2f 



< z<w. 



If the thickness w is small enough so that the last term can be neglected, we obtain, 



2 \2 



nXf 



n 2 X 
16 



(30) 



(31) 



which is the conventional zone plate law for 1:1 imaging^. The inset on bottom of Fig. |l2| shows the zone plate 
structure that satisfies Eq. (31 ) (solid flat lines), and the ideal one according to Eq. (30) (dashed curved lines) for 1:1 
imaging. In this case, the flat zones show an increased deviation from the ideal curved ones with increasing radius, 
indicating there may exist a theoretical limit preventing the NA from approaching unity even for 1:1 imaging. To 
investigate this limit, we studied the imaging property of a flat MLL when illuminated by a spherical wave emitted 
from a point source 2/ away from its center (Fig. 12 1. The parameters for the MLL are: w = 13.5 /im, / = 0.472 



mm, x max = 75 /im and the outmost zone width of 0.2 nm. It is easy to derive the phase function for this MLL 
structure from Eq. (311, 



<f> h = Ahk(y/x 2 /A + f 2 -f). 



By solving Eq. ( 12 ), we are able to simulate the isophotes near the image point (Fig. 12 1. We found a peak FWHM of 



0.53 nm, 33% larger than the diffraction limit from Rayleigh criterion, 0.4 nm. The larger focused beam size, compared 
to the diffraction limit, is a result of deviating from the ideal structure. From the drawing of MLL structures in Fig. 
[L2| as the radius increases the flat zone profile deviates further from the ideal elliptical shape. Therefore, as opposed 
to the conclusion made by Pfeiffer et alP^l, we argue that according to our calculation there is a resolution limit even 
for 1:1 imaging when an MLL with flat zone profiles is used. 



V. SUMMARIES AND CONCLUSIONS 



In summary, we present a formalism of dynamical diffraction from Fresnel zone plates analogous to the Takagi- 
Taupin equations for a strained single crystal. The basic equations for dynamical diffraction are derived and employed 
to study diffraction properties of Multilayer-Laue-Lense with various types of zone profiles, including flat, tilted, 
wedged and ideally curved. Our study shows that for a thick MLL made of flat zones, in case of plane wave illumination 
with normal inclination angle, the dynamical diffraction effect of the zone plate would prevent it from focusing the 
beam size down to a few nanometers. A rough estimation in this case shows that when the thickness is optimized for 
dynamical diffraction, the achievable NA is about V2An 7 a result similar to that for waveguides, and in an agreement 
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with the conclusion obtained recently by Pfeiffer et afl^. We also demonstrate that the achievable NA of an MLL 
with flat zones can be increased by a factor of roughly 2 by tilting zones. The MLL with wedged zones, which satisfy 
locally the Bragg condition, is studied as well. It is observed that wedged MLL's can reach a NA of ~0.1 when high 
diffraction efficiency is achieved, i.e., effectively focusing hard x-rays down to several angstroms. Ideal structures 
which satisfy both the Bragg condition and the phase requirement are obtained with the aid of an Ewald sphere 
construction. Not unexpectedly, they turn out to be confocal paraboloids for an incident plane wave and confocal 
ellipsoids for an incident spherical wave. 
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Flat Titled Wedged Curved 




FIG. 1: Types of multilayer-Laue-lens (MLL) layer arrangements: flat, tilted, wedged, and curved. 
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FIG. 2: A mapping of layers between a Fresnel zone plate and a periodic structure composed by binary bars with thickness 
ratio 1:1. 




FIG. 3: The zone plate structure can be represented by a set of sine waves with varying period that obey zone plates law. 
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FIG. 5: a) The radial efficiency distributions of different diffraction orders on the exit surface of a fiat MLL at 6 — (top), 
8 = 1.6 (middle) and — 3.2 mradian (bottom), b) The efficiency variation of the —I s * order diffraction inside the zone plate 
at 8 = (left) and 8 = 1.6 mradian (right), c) The radial efficiency distribution of the —I s order diffraction on the exit surface 
at different thickness and tilting angles. 



1G 




FIG. 6: a) Isophotes near the focus at different tilting angles. 8 = 0, 1.1, 2.1 and 3.2 mrad. in i), ii), iii) and iv). The intensity 
is plotted in logarithm scale, b) The intensity profiles on the best focal plane. 
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FIG. 7: a) The radial intensity distribution of the -1st order diffraction on the exit surface of a wedged MLL with an outmost 
zone width of 0.25 nm. Two-beam approximation is assumed. The dashed curve represents the phase deviation from a perfect 
converging spherical wave front, b) The intensity profile on the focal plane, showing a peak width FWHM=0.34 nm. The inset 
on top is the isophotes near focus, intensity in logarithm scale. The one on bottom is a sketch of the wedged MLL structure. 
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FIG. 9: a) A schematic of satisfying Braggs law everywhere for an incident plane wave. All diffracted wave vectors point to 
the focal point F. b) A structure consisting of a set of confocal paraboloids obtained from a), serving as an ideal structure to 
focus an plane wave into a point. 
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FIG. 10: The radial intensity distribution of the —I s * order diffraction on the exit surface of an MLL with ideal parabolic zones 
and outmost zone width 0.25 nm, based on the 2-beam simulation. The dashed curve represents the phase deviation from a 
perfect converging spherical wave front, b) The intensity profile on the focal plane, showing a peak with FWHM=0.21 nm. 
The inset on top is the isophote pattern near focus. Intensity is on a logarithm scale, and the one on bottom is a sketch of the 
ideal structure 
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FIG. 11: a) A schematic of satisfying Bragg's law everywhere for an incident spherical wave. All diffracted wave vectors point 
to the focal point F. b) A structure consisting of a set of confocal ellipsoids, serving as an ideal structure to focus a point in 
object space into a point in image space. 
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FIG. 12: The focused beam profile on the plane of best focus and the isophotes near the image point (the insert on top) for 
1:1 imaging. The insert on bottom shows the difference of an MLL with fiat zones (solid line) and an ideal MLL (dashed line). 



